Decoupling phenomena in supercooled liquids: Signatures in the energy landscape 
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A significant deviation from the Debye model of rotational diffusion in the dynamics of orienta- 
tional degrees of freedom in an equimolar mixture of ellipsoids of revolution and spheres is found 
to begin precisely at a temperature at which the average inherent structure energy of the system 
starts falling with drop in temperature. We argue that this onset temperature corresponds to the 
emergence of the a-process as a distinct mode of orientational relaxation. Equally important, we 
find that the coupling between the rotational and translational diffusion breaks down at a still lower 
temperature where a sharp change occurs in the temperature dependence of the average inherent 
structure energy. 
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The relaxation phenomena in supercooled liquids con- 
tinue to stimulate intense research interests despite per- 
sistent research activity over decades 0, 0- A major 
contribution to this activity has come from a variety of 
experimental techniques, e.g., dielectric relaxation spec- 
troscopy 0, light scattering 0, time resolved optical 
spectroscopy NMR spectroscopy 0, 0, and opti- 
cal Kerr effect spectroscopy ||, that probe dynamics 
of orientational degrees of freedom (ODOF). These ex- 
periments reveal an array of dynamical features some of 
which might owe their origin to the nontrivial interplay 
between orientational and translational degrees of free- 
dom On the other hand, a decoupling between ro- 
tational and translational diffusion is observed in deeply 
supercooled molecular liquids in the sense that orienta- 
tional correlation time continues to track the viscosity as 
given by the Debye-Stokes-Einstein (DSE) relationship 
while translational diffusion coefficient does not in con- 
trary to what is predicted by the Debye-Einstein (DE) 
relation [H [ll|. The a - (3 bifurcation d E 0, 
which commonly refers to the bifurcation of two peaks 
in the dielectric relaxation spectra marks yet an- 
other decoupling, this time between two distinct mech- 
anisms for orientational relaxation in liquids composed 
of non-spherical molecules. The bifurcation temperature 
Tb was believed to be close to the mode-coupling criti- 
cal temperature T c but it has been recently shown 
in broadband dielectric relaxation measurements that a- 
and /3-relaxations merge together only well above T c • 

Here we address the decoupling phenomena from the 
perspective of potential energy landscape by studying a 
system with orientational degrees of freedom across the 
supercooled regime. The energy landscape formalism is 
an approach that explores the features of the underlying 
potential energy surface of a system for understanding its 
complex dynamics [la, UJ1 lla lii! Ui| • In an appealing 
landscape study of a binary mixture of spheres [2fJ , Sas- 
try et. al. found that the onset of non-exponential relax- 
ation in the supercooled regime corresponded to the tem- 
perature, below which the dynamics of the system was 



influenced by its energy landscape. As in Ref. |20j, the 
focus in most molecular dynamics simulation studies on 
supercooled liquids, with a few notable exceptions [2l| . 
has been on atomic systems which involve only transla- 
tional degrees of freedom (TDOF) H3- 

We investigate a binary mixture of 128 Gay-Berne 
ellipsoids of revolution and 128 Lennard- Jones spheres 
along an isochore at a series of temperatures down to 
the deeply supercooled state |2^|. The choice of such a 
system is motivated by the success of binary mixtures of 
Lennard- Jones spheres |2^] , that are widely used for com- 
puter simulations of supercooled liquids. In our system 
|22j , the interaction between the ellipsoids of revolution is 
given by the Gay-Berne pair potential |25j , which explic- 
itly incorporates anisotropy in both the attractive and 
the repulsive parts of the interaction with a single-site 
representation for each ellipsoid of revolution. The in- 
teraction potential between a sphere and an ellipsoid of 
revolution, which is chosen to be a prolate (with aspect 
ratio k = 2), is given by following Cleaver and coworkers 
|26j . We have determined the set of energy and length 
parameters such that neither any phase separation occurs 
nor any liquid crystalline phase with orientational order 
appears even at the lowest temperature studied at a high 
density HjJI^J. Across the supercooled regime, the dif- 
fusion coefficients for both the spheres and the ellipsoids 
of revolution are found to follow a power law temperature 
dependence, i.e., D t = C D (T - T C )^ D , with T c = 0.454 
for the former and T c — 0.460 for the latter. This sug- 
gests, within the error limit, the critical temperature T c 
to be independent of the type of particle as predicted by 
the mode-coupling theory (MCT). 

In Fig. 1, we show the temperature dependence of the 
average energy of the inherent structures for our binary 
system. At high temperatures (T > 1.0), the average in- 
herent structure energy remains fairly insensitive to tem- 
perature variation. Below T ~ 1.0, this energy decreases 
progressively up to the lowest temperature studied here. 
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FIG. 1: The average inherent structure energy per particle 
eis of our binary system as a function of temperature over 
the whole temperature range studied. The inset plots eis 
versus 1/T over the temperature range across which the aver- 
age inherent structure energy is on a decline. The solid lines 
are the linear fits to the data over two distinct temperature 
regimes. 



We find that this crossover temperature corresponds to 
the onset of non-exponential relaxation in the decay of 
the self intermediate scattering function F s (k,t) (data 
not shown). This can be taken as a dynamical signa- 
ture of the crossover behavior in translational degrees of 
freedom. Such a crossover in the average depth of the po- 
tential energy minima explored by the system has been 
observed earlier , and shown to be consistent with the 
thermal sampling of a Gaussian distribution of energies 
for the local minima 19, 28] . The latter, within harmonic 
approximation, predicts an inverse temperature depen- 
dence of the average inherent structure energy [28|. We, 
however, observe two distinct temperature regimes, where 
the inverse temperature dependence of the average inher- 
ent structure energy holds true, with a sharp change at 
a second crossover temperature T ~ 0.6 as illustrated in 
the inset of Fig. 1. We now investigate whether the tem- 
perature dependence of the average inherent structure 
energy has any correlation with change in the dynamics 
of orientational degrees of freedom. 

Figure 2a shows how the ratio of the first to second 
rank rotational correlation times, t±/t2, evolves as tem- 
perature drops. The Z-th rank rotational correlation time 
t; is defined as t; = C^^dt, where c\ s \t) is the l- 
th rank single-particle orientational correlation function: 
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FIG. 2: a. The temperature dependence of n(T)/n(T), the 
ratio of the first to second rank rotational correlation times 
(squares). On a different scale (appearing on the right) shown 
again is the temperature dependence of the average inherent 
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ison (circles), b. The time evolution of the single-particle 
second-rank orientational correlation function at several tem- 
peratures with a shift in the time origin to to = 1.2 followed 
by normalization to the corresponding value at to- C. The ex- 
ponents values /3i(T) and 02 (T), obtained from the fits to the 
stretched exponential form of the first and second rank single- 
particle orientational correlation functions, respectively, as a 
function of temperature. 
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FIG. 3: The product of the translational diffusion coefficient 
Df and the second-rank rotational correlation time T2 for 
the ellipsoids of revolution as a function of temperature. In- 
set: The inverse temperature dependence the logarithm of 
the second-rank rotational correlation time T2 (circles). On 
a different scale (appearing on the right of the inset) shown 
is the inverse temperature dependence the logarithm of the 
shear viscosity (squares). The solid and dashed lines are the 
respective Arrhenius fits to data over a restricted temperature 
range. 



metry of the ellipsoid of revolution i, Pi is the Ith rank 
Legendre polynomial and the angular brackets stand for 
ensemble averaging. It is evident in Fig. 2a that the ratio 
has a value close to 3 as predicted by the Debye model of 
rotational diffusion at high temperatures (T > 1.0) and 
starts declining steadily from T ~ 1.0 until it reaches 
a value nearly unity at low temperatures. In the diffu- 
sive limit, the orientational motion occurs in small steps, 
while a value for this ratio close to 1 is taken to suggest 
that the orientational motion takes place via long angular 
jumps [l^ . 

At high temperatures, the long time decay of (t) 
is exponential while the Kohlrausch- Williams- Watts 
stretched exponential form provides a reasonable fit to 
the long-time behavior at low temperatures. To elimi- 
nate the short-time Gaussian time dependence, we show 
in Fig. 2b the time evolution of the function C% (t — 

to)/C'2 S \to) for t > to, and consider the stretched expo- 
nential form exp[-((t - £ )/t(T))' 32(t) ] that takes into 
account this transformation to fit the data. The devia- 
tion of the exponent /3 2 (T) (0 < (3 2 {T) < \) from unity 
is a measure of the degree of non-exponential relaxation. 
Figure 2c shows the temperature dependence /^(T) and 
also (3i(T), the latter corresponding to the long-time de- 
cay of C[ s) (t) (not shown here). While both /? X (T) and 
foiT) are very close to unity at high temperatures, they 



FIG. 4: The translational diffusion coefficient Dt versus the 
shear viscosity 77 in a log-log plot for both spheres (circles) and 
the ellipsoids of revolution (squares). The solid and dashed 
lines are the respective linear fits. 

start dropping as temperature falls. It is evident that the 
stretching is more pronounced in the second rank orien- 
tational correlation function as compared to the first rank 
orientational correlation function. Such an observation 
has been reported previously while comparing data from 
dielectric and NMR experiments [jj. We note that the 
signature of non-exponential relaxation in C% (t) first be- 
comes appreciable and later gets progressively more pro- 
nounced as temperature drops below the onset tempera- 
ture T ~ 1.0. 

It follows from above that the onset of the growth of 
the depth of the potential energy minima explored by the 
system correlates with a change in the mechanism of ori- 
entational motion from being simply diffusive. There is 
evidence for the /^-relaxation to be diffusive in character 
and the orientational relaxation above Tb is associated 
with the /3-process only 0. Thus, the onset tempera- 
ture corresponds to the emergence of the a-process as a 
distinct mode of orientational relaxation and appears to 
be coinciding with the bifurcation temperature Tb- We, 
however, find that the latter is somewhat higher than T c 
as indeed observed in Ref. |15| in contrary to what often 
believed. Nevertheless, Fig. 2a illustrates that the decline 
of the ratio Txjr-i closely tracks the fall of the average in- 
herent structure energy of the system. Such a correlation, 
to the best of our knowledge, has not been demonstrated 
before. Stillinger interpreted the a — [3 bifurcation in 
terms of the topography of the potential energy land- 
scape 0. In Stillinger's picture, the /3-processes cor- 
respond to the elementary relaxations between contigu- 
ous basins while the a-processes invoke escape from one 
metabasin and eventually into other with an involvement 
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of high free energy of activation. Such a description is 
consistent with the growth of the depth of the potential 
energy minima explored by the system below the bifur- 
cation temperature. 

We now address the decoupling between rotational and 
translational diffusion. The combination of the SE and 
DSE equations predicts the product D t T2 to be indepen- 
dent of temperature even when the macroscopic observ- 
able viscosity increases by many orders of magnitude as 
the glass transition temperature T g is approached from 
above [10( . Figure 3 shows that such a relationship breaks 
down at T ~ 0.6 and below with the product growing fast 
with decrease in temperature. The inset of Fig. 3 illus- 
trates that the decoupling between the two microscopic 
observables occurs precisely at the same temperature at 
which both the orientational correlation time r 2 and the 
viscosity start showing steady deviation from the Arrhe- 
nius temperature behavior. The inset of Fig. 1 shows 
that at this temperature the linear variation of the av- 
erage inherent structure energy with the inverse temper- 
ature undergoes a sharp change with an increase in the 
rate of fall. 

In Fig. 4, we show the variation of the translational dif- 
fusion coefficient D t with the coefficient of shear viscos- 
ity i] in a log-log plot over the whole temperature range 
studied here. The linearity of the curve implies a power 
law dependence: D t oc r)~ a , a being the exponent, for 
both the spheres and the ellipsoids of revolution. We 
find a = 0.83 for the former and a = 0.75 for the lat- 
ter. The fractional power law dependence suggests the 
enhancement of translational diffusion relative to what 
the SE relationship predicts. The a values obtained here 
compare well with 0.77, observed by Ediger and cowork- 
ers in a direct measurement of self-diffusion of a single- 
component glass- forming liquid reported recently 29] . 

We now argue that the deviation from the Debye be- 
havior of rotational diffusion is due to a crossover from a 
collision dominated to a correlation dominated regime of 
orientational relaxation. Below T ~ 1, the equilibrium 
pair correlations are expected to grow as e/s starts its 
descend. The effects of correlations can be included via 
the standard Zwanzig-Mori continued fraction represen- 
tation: Ci{z) = l/(z + l(l + l)k B T/(I(z + fi{z))), where 
the memory kernel F;(z) can be decomposed in the spirit 
of MCT as ti(z) ~ f coU {z) + AT t (z), the effect of corre- 
lations being contained in ATi(z), and z stands for the 
Laplace frequency. The rank dependence of the memory 
function can be approximately expressed in terms of the 
torques-torque time correlation function |3fJ. Because 
of the up-down symmetry of the ellipsoids of revolution, 
the contribution of the correlation to Ar 2 is expected to 
be larger than that to Afi, particularly at t = 0. Thus, 



f 2 > f i and this would lower the value of the ratio ri/T2. 
An opposite trend has been observed for dipolar systems 
where l\ > T 2 and an upward deviation for the ratio 
Ti/r 2 is observed [3fil |. 

In summary, we have established a correlation of the 
breakdown of the Debye model of orientational relax- 
ation and the a — [3 bifurcation with the manner of ex- 
ploration of the underlying potential energy landscape 
in a model system. Equally important, the decoupling 
between the rotational and translational diffusion is sig- 
naled by a sharp rise in the rate of fall of the average 
inherent structure energy with the inverse temperature. 
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